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Forward 

More  than  two  decades  ago  Bahar  formulated  a  full  wave  theory  for  rough  surface 
scattering(l,2).  This  theory  was  developed  along  the  lines  used  for  tapered  waveguides 
and  based  on  the  use  of  generalized  telegraphists  equations.  As  a  first  step  Bahar 
introduced  a  spectral  expansion  of  the  fields  in  the  transverse  plane  and  in  the  direction 
perpendicular  to  the  interface  separating  two  different  media.  This  formulation  is 
essentially  limited  to  treating  rough  surfaces  that  are  rough  in  one  direction  only,  although 
Bahar  did  extend  the  theory  in  a  more  or  less  ad  hoc  way  to  treat  rough  surfaces  with  two 
dimensional  roughness  governed  by  isotropic  statistics  (3).  By  expanding  the  fields  in  the 
vertical  direction  in  terms  of  the  spectral  expansion  and  substituting  into  Maxwell’s 
equations  one  obtains  certain  coupling  integrals  that  corresponded  to  mode  coupling  in  the 
analogous  waveguide  problem,  along  with  some  boundary  terms.  The  resultant  equations 
are  too  difficult  to  solve  exactly  so  Bahar  resorted  to  an  iterative  method.  In  order  to 
obtain  a  first  iterative  solution  Bahar  approximates  the  fields  in  the  mode  coupling 
integrals  by  an  approximate  solution  corresponding  to  a  plane  wave  incident  and  reflected 
from  the  surface  using  the  approximation  that  locally  the  field  is  reflected  as  though  the 
rough  surface  was  flat  with  a  height  equal  to  the  local  height.  Later  on  Bahar  modified  his 
theory  by  referring  the  fields  and  scattering  coefficients  to  the  local  normal.  A  detailed 
study  of  this  theory  has  shown  that  it  has  serious  limitations  and  is  generally  unreliable 
despite  claims  made  by  Bahar  as  to  its  wide  applicability  (4,5).  Since  the  first  order  field 
that  Bahar  uses  for  the  first  iterative  solution  is  in  itself  not  a  solution  to  Maxwell’s 
equations  Bahar’ s  method  is  not  a  systematic  iterative  or  perturbation  solution.  For  this 
reason  it  is  difficult  to  proceed  to  higher  order  solutions  since  the  starting  point  already 
has  errors  that  are  difficult  to  quantify. 

The  above  somewhat  shaky  foundation  for  the  full  wave  theory  led  us  to 
reconsider  the  basic  theory  from  the  begiiming  and  to  follow  a  more  rigorous 
mathematical  derivation  as  far  as  it  could  be  pursued  without  approximations.  The  first 
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step  was  the  recognition  that  a  spectral  expansion  of  the  field  components  in  the  vertical 
direction  above  the  rough  surface  could  be  carried  out  using  basis  functions  that  locally 
satisfy  the  boundary  conditions  for  a  surface  that  is  flat  but  having  the  local  height  (6).  By 
using  this  expansion  in  Maxwell’s  equations  one  can  formulate  a  full  wave  theory  for 
rough  surfaces  that  have  arbitrary  roughness  in  two  directions,  without  any  need  to 
introduce  the  telegraphists  equations.  We  have  also  found  that  it  is  possible  to  carry  the 
theory  forward  without  approximations  to  a  point  where  no  mode  coupling  integrals  are 
present,  the  only  unknowns  being  boundary  values  for  the  unknown  field.  We  call  this 
method  the  local  spectral  expansion  method  to  distinguish  it  from  Bahar’s  method  which  is 
based  on  the  telegraphists  equations.  Based  on  this  approach  it  was  of  interest  to  examine 
what  results  could  be  obtained  if  the  same  tangent  field  approximation  to  the  surface  fields 
as  used  in  the  Kirchhoff  method  was  also  used  in  the  local  spectral  expansion  method.  We 
found  that  this  led  to  useful  results  that  had  a  range  of  application  that  covered  both  the 
small  perturbation  regime  as  well  as  the  Kirchhoff  regime.  A  curious  aspect  of  this  more 
rigorous  method  was  that  the  solutions  obtained  were  intrinsically  non-reciprocal,  similar 
with  what  happens  in  the  phase  perturbation  method  (7).  In  spite  of  this  lack  of  reciprocity 
the  predicted  scattering  cross  sections  turn  out  to  be  more  accurate  than  those  found  using 
either  the  small  perturbation  solution,  the  Kirchhoff  solution,  or  Bahar’s  full  wave 
solution.  Although  reciprocity  is  physically  required, reciprocity  or  the  lack  of  it  in 
approximate  solutions  is  not  a  criteria  for  the  accuracy  of  the  approximate  solution. 

A.  Statement  of  the  Problem  Studied. 

The  problem  investigated  was  the  use  of  the  local  spectral  expansion  method  along 
with  the  Kirchhoff  tangent  field  approximation  to  the  surface  fields  to  evaluate  the 
scattering  cross  sections  of  rough  surfaces  governed  by  Gaussian  statistics,  and  to 
compare  the  results  with  numerical  solutions  based  on  the  Monte  Carlo  method  as  well  as 
those  predicted  by  other  theories.  Perfectly  conducting  surfaces  with  roughness  along  one 
direction  as  well  as  along  two  directions  were  studied.  In  addition  scattering  from  rough 
dielectric  surfaces  were  also  investigated,  but  limited  to  surfaces  with  one  dimensional 
roughness.  This  work  was  to  be  used  to  develop  higher  order  iterative  solutions.  The 
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foundations  for  higher  order  iterative  solutions  were  developed  in  a  mathematically  sound 
manner  but  no  results  beyond  using  the  first  order  Kirchhoff  tangent  field  approximation 
was  obtained  during  the  course  of  the  present  investigation. 


B.  Summary  of  the  Most  Important  Results. 

Consider  a  perfectly  conducting  rough  surface  with  one  dimensional  roughness 
described  by  the  height  profile  h(x).  The  y  axis  is  perpendicular  to  the  mean  surface  which 
coincides  with  the  y=0  plane.  For  TE  and  TM  waves  the  following  local  basis  functions 
are  used 


V«(x.y.ky) 


It  is  convenient  to  define  the  following  transform  of  the  scattered  electric  field  for  TE 
waves  00 


A»(x,ky)H-^  fE*(x.y)e^»My 


The  transform  of  the  scattered  electric  field  in  terms  of  the  odd  local  basis  function  can  be 
expressed  in  terms  of  the  above  transform.  The  equation  satisfied  by  the  A®  transform  was 
found  to  be 


a*  A* 

ax» 


+(k»-k;)A*  = 


_1 _ ^ 
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Note  that  the  only  unknown  in  this  equation  is  the  normal  derivative  of  the  scattered 
electric  field.  This  normal  derivative  is  related  to  the  unknown  current  density  on  the 
surface,  i.e.  the  tangential  magnetic  field  at  the  surface.  In  the  tangent  field  approximation 
this  unknown  normal  derivative  is  approximated  by  the  normal  derivative  of  the  incident 
plane  wave  electric  field.  The  resultant  equation  can  now  be  formally  solved  using  a 
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Fourier  transform  with  respect  to  the  x  coordinate.  From  this  solution  the  scattered 
electric  field  in  the  far  zone  can  be  obtained  by  using  the  method  of  stationary  phase. 
This  solution  is 


E®(x,y)  =  -^^^j=^=^Fi(0i,0s)  J -  l)dx' 


-L 

L 


-jk.  r+JT 


-L 

L 


where  li,=  a,k* +tyk*  =  k,(a,sin0,+ayCOS0^,  v,  =  k',  -  k*,  =  k^(sinei-sin0J, 

Py  =  k‘y  +  k;  =  k„(cos0j+cos0J,  qy  =  k‘y  -  k*y  =  k^CcosGi-cose^,  Fi  =  (v J  +  pj  ) /py  = 

2ko(l  +  cos(0i  +e,))/(cos(ei)  + 008(0, )),  F2  =(vj  +<ly)/<ly  “ 

2ko  (1  -  cos(0i  -  08))/(cos(0i )  -  008(0,)) 

The  scattering  cross  sections  can  be  evaluated  using  standard  procedures  when  it  is 
assumed  that  the  surface  roughness  is  governed  by  Gaussian  statistics. 

Some  typical  results  that  have  been  obtained  are  shown  in  Fig.’s  1-3  and  are 
compared  with  results  obtained  using  the  small  perturbation  method,  the  Kirchhoff 
method,  Bahar’s  regular  full  wave  method,  and  numerical  solutions  provided  to  us  by  Dr. 
E.LThorsos.  The  results  obtained  with  the  local  spectral  expansion  method  (LSEM)  and 
shown  in  Fig.’s  1  and  2  show  good  agreement  with  the  numerical  evaluations  and  are 
clearly  superior  to  those  obtained  using  the  regular  full  wave  method  (RFWM)  or  the 
Kirchhoff  method  (KA).  In  Fig.  3  the  angle  of  incidence  is  80°  and  for  this  near  grazing 
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angle  of  incidence  none  of  the  approximate  theories  are  accurate  in  the  back  scatter 
direction. 

Similar  results  were  also  obtained  for  TM  polarization  but  there  was  less  in  the 
way  of  numerical  solutions  with  which  to  compare  the  results.  For  TM  polarization  the 
local  spectral  expansion  method  appears  to  be  closer  to  the  Kirchhoff  results. 

Computations  of  the  scattering  cross  sections  for  rough  dielectric  surfaces  were 
also  carried  out  and  compared  with  the  numerical  results  of  Sanchez-GU  and  Nieto- 
Vesperinas  (8)  for  the  case  of  an  angle  of  incidence  of  40°  and  a  relative  permittivity  of 
2.04.  The  agreement  with  the  numerical  simulations  are  better  than  those  obtained  using 
the  Kirchhoff  approximation  or  the  regular  full  wave  theory,  as  can  be  seen  from  Fig.  4. 
Again  there  is  a  scarcity  of  numerical  results  to  compare  with. 

Overall  it  was  found  that  by  using  the  tangent  field  approximation  with  the  local 
spectral  expansion  method  a  more  accurate  prediction  of  the  scattering  cross  sections 
could  be  made  than  that  provided  by  either  the  small  perturbation  method  or  the  Kirchhoff 
methods  alone.  It  is  generally  believed  by  the  current  investigators  that  the  local  spectral 
expansion  method  is  a  useful  contribution  to  the  collection  of  approximate  theories  for 
rough  surface  scattering. 
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Scattering  angle  (deg.) 

Rg.  1,  Scattering  cross  section  for  horizontal  polarization  (TE  waves)  for  a  perfectly 
conducting  rough  surface.  Angle  of  incidence  =  45®,  mas  surface  height  0.212  wavelength, 
correlation  length  0.9  wavelength,  rms  slope=  18.42°. 


Scattering  angjle  (deg.) 

Hg.  2.  Scattering  cross  section  for  horizontal  polarization  (TE  waves)  for  a  perfectly 
conducting  rough  surface.  Angle  of  incidence  =  60®,  rms  surface  height  0.048  wavele 
correlation  length  0.796  wavelength,  rms  slope=4.84®. 


Scattering  angle  (deg.) 

Fig.  3.  Scattering  cross  section  for  horizontal  polarization  (TE  waves)  for  a  perfectly 
conducting  rough  surface.  Angle  of  incidence  =  80®,  rms  surface  height  0.107  wavelength, 
correlation  length  0.45  wavelength,  rms  slope=18.59®. 


Normalized  diffuse  intensity 


Scattering  angle  (deg.) 


Fig.  4.  Scattering  cross  section  for  a  rough  dielectric  surface.  Angle  of  incidence  =  40°, 
relative  permittivity  =  2.04,  rms  surface  height  0.5  wavelength,  correlation  length  3.16 
wavelength. 


